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Abstract 

a: 

Ph! The transverse polarization distribution of quarks hi(x) is computed in a confine- 

ment model, the chiral chromodielectric model. The flavor structure of hi, its Q 2 
evolution and Soffer's inequality are studied. The Drell-Yan double transverse asym- 
metry Att is evaluated and found to be one order of magnitude smaller than the double 
longitudinal asymmetry. 
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The interest in the transverse polarization distribution of quarks and antiquarks, 
customarily called hi (x) , has been recently strengthened by the perspective of its pos- 
sible measurement in future collider experiments. Originally introduced by Ralston 
and Soper [|J, who called it /it(^), h\(x) has been studied in detail, from a formal 
point of view, in more recent papers 0, [|. The possible ways of measuring h\(x) in 
various hard processes have been also thouroughly investigated f|, ^ |6|]. Despite 
this intensive work, not much is actually known about the shape and the magnitude 
of hi(x). This is of course not surprising since hi, like all quark distributions, cannot 
be derived from the fundamental theory of strong interactions, QCD. At present we 
possess only: i) an admittedly crude evaluation of hi(x) in the simplest version of the 
MIT bag model ||, and ii) an estimate of its first moment - the so-called tensor charge 
- obtained by QCD sum rule methods 0. It is clear that a more sophisticated model 
calculation of hi(x) is called for. This would also provide useful indications about the 
concrete possibility of a measurement of hi(x) in the experiments which are now being 
planned. 

The difficulty of an experimental determination of hx(x), which is a leading twist 
quantity, resides mainly in the fact that, being a chirally-odd distribution, it is not 
measurable in polarized deep inelastic scattering. The best method to extract hi{x) 
seems to be j|, [IJ the Drell-Yan dilepton production with two transversely polarized 
proton beams, an experiment which will be performed in the near future at RHIC M . 
The Drell-Yan double transverse asymmetry Att contains information on the flavor 
structure of h\. Therefore it would be important to predict the magnitude of Att i n 
order to test the feasibility of the experiment. The available model calculations of Att 
|| g] are all based on the assumption that h\ is approximately equal to the helicity 
distribution Aq. Although this is probably true at very low momentum scales, such 
as those at which confinement model computations are implicitly performed (Q 2 < 1 
GeV 2 , see below), it is certainly not true at experimental Q 2 scales (say Q 2 > 10 GeV 2 ). 
The reason is that Aq and h\ evolve differently in Q 2 . Whereas the first moment of Aq 
is constant, the first moment of h\ decreases with increasing Q 2 , and the evolution in 
the x-shape is even more dramatically different. Again, a more firmly based evaluation 
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of Att is needed to check whether this quantity is comparable in magnitude with the 
double longitudinal asymmetry A LL , as it is claimed in || [|. 

Another issue which is certainly worth exploring in the framework of confinement 
models is the inequality among leading twist polarized and unpolarized distribution 
functions recently derived by Soffer (see also [[U]]) and the possibility of its satu- 
ration. 

In the following we shall provide a theoretical determination of hi (x) in a confine- 
ment model, the chiral chromodielectric model |I2|, which has been already succesfully 



used to compute other leading twist structure functions and various nucleon properties. 
In particular, the flavor structure of hi will be described in detail and the effects of the 
peculiar QCD evolution of hi investigated. A prediction for A TT will also be presented. 
As we shall see, due to the different evolution of h\ and Aq, Att turns out to be much 
smaller than A LL . 

The quark transverse polarization distribution reads || 

hi(x) = / dre"^ +r (NS ± \^i(Ol±lsMO)\NS±)k + =t ± =o. (1) 
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A similar expression holds for the antiquark distribution, with the exchange of ip and ip^ ■ 
Note that in eq. ([J) only the 'good' light-cone components of the fields, ip+ — 
appear, signaling that hi is a leading twist quantity. The hi distribution measures the 
difference in the number of quarks with transverse polarization parallel (|) and anti- 
parallel (|) to the proton transverse polarization. This can be made transparent by 
introducing the Pauli-Lubanski projectors P]_ = |(1±7j_75) and inserting a complete 
set of states {|A)} in eq. (p. We then get 

Mo;) = 4E{KiVSjPi^0)^ (2) 
V '2 x 

In a (projected) mean-field approximation, the matrix elements in eq. @ can be 
rewritten in terms of single-particle (quark or antiquark) matrix elements. For a flavor 
/ one thus gets 

* w = 75 E E nf, a, m) / amw-W- pi] 

x Jp(p a , m)7 + 7 ± 7 5 y?(p a , m) , (3) 



where (p is the single-quark wave function, m is the projection of the quark spin along 
the direction of the nucleon's spin, P(f,a,m) is the probability of extracting a quark 
of flavor / and spin m leaving a state generically labelled by the quantum number a. 
The overlap function A a (p a ) contains the details of the intermediate states and of the 
projection used to obtain a nucleon with definite linear momentum from a three-quark 
bag (see for instance |L3], |T^|). The intermediate states which contribute to eqs. 
are 2q and 3qlq states for the quark distribution, and 4q states for the antiquark 
distribution. 

At this point, we can already discuss qualitatively the Soffer inequality |]I0 |, which 
reads 

q f (x) + Aq f (x) > 2\h({x)\ , (4) 

where Aqf is the helicity distribution function and qf the unpolarized density. This 
relation has been proved in the parton model [|TI], |TT|] (for a QCD-improved parton 
model discussion of Soffer's inequality see |TB[), and is satisfied flavor by flavor by both 
the quark and the antiquark distributions. An interesting issue is whether Soffer's 
inequality is saturated in some quark model (which means that \hi | takes its maximal 
value). To clarify this problem let us write the various leading twist distributions in 
an explicit form 



a m 

= EE p (/>«. m )(-i) (m+3/2) ^) 

a m 

h{{x) = EE p (/.«.™)(-i) (m+3/2) ^W 



(5) 
(6) 
(7) 



where 




dp* 



(2vr)3(2p0 



-A Q (p a )5[(l - x)p + -p+] 



1 

x - < 
2 



u 2 (p a ) + 2u(p a )v(p a )^- l + v 2 (p a ) 



u\p a ) + 2u(p a )v(p a ) J ^- l + v\p a ){2 (j^y 
u 2 (p a ) + 2u(p a )v(p a )^ + v 2 {p a )(l - (fifr 
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In eq. (§) p is the projection of the momentum in the plane perpendicular to the 
proton's trajectory (chosen to be the z axis), and the currents have been written in 
terms of the single quark wave-function in momentum space 

^ m > = < 9 > 

Notice that the three quantities F a ,G a , H a satisfy the equality: F a (x) + G a (x) = 
2H a (x). This has led to the erroneous conclusion |L0| that the inequality is sat- 
urated for a relativistic quark model, such as the MIT bag model. It is clear from 
eqs. (5-7) that the spin-flavor structure of the proton, which results in the appear- 
ance of the probabilities P(f,a,m), spoils this argument and prevents in general the 
saturation of the inequality 

Soffer's inequality is saturated only in very specific (and somehow unrealistic) cases. 
For instance, it is saturated when P(f,a, —1/2) = 0, which happens if the proton is 
modeled as a bound state of a scalar diquark and a u quark. Of course, this is too a 
rough picture of the proton. However, it is interesting to note that in 577(6) the A is 
a bound state of a scalar-isoscalar ud diquark and an s quark: the h\ distribution of 
the latter then attains the maximal value compatible with (|j). 

Another instance of saturation is when F a = G a = H a and P(f, a, — 1/2) = 
2 P(f, a, 1/2). It is easy to verify that this happens for the d quark distribution in 
a nonrelativistic model of the proton with an SU (6) wavefunction. 

Apart from the two particular cases illustrated above, Soffer's inequality should not 
be expected to be saturated, and indeed it is satisfied but not saturated in the model 
we present here. 

The model of the nucleon that we use to compute all ingredients appearing in eq. 
and then to evaluate the quark distribution functions is the chiral chromodielectric 
model (CCDM) 0. The Lagrangian of the CCDM reads 

C = i^^dai) + —1p (<J + Z75T • 7r) "0 

X 

+ ^(^x) 2 -^MV + ^(^) 2 + i(^7r) 2 -f/(a,7r), (10) 

where U(a, 7r) is the usual mexican-hat potential, see e.g. fll6"| . C describes a system 
of interacting quarks, pions, sigmas and a scalar-isoscalar chiral singlet field x- The 
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parameters of the model are: the chiral meson masses = 0.14 GeV, m a = 1.2 GeV, 
the pion decay constant f n = 93 MeV, the quark-meson coupling constant g, and the 
mass M of the x field. The parameters g and M, which are the only free parameters of 
the model, have been univoquely fixed by reproducing the average nucleon-delta mass 
and the isoscalar radius of the proton. 

The CCDM Lagrangian ( |I~0"D contains a single-minimum potential for the chro- 
modielectric field x- V(x) = \M 2 x 2 - A double-minimum version of the CCDM is also 
widely studied and used (see for instance [17]). We have checked that the structure 
functions computed in the two versions of the CCDM do not differ sensibly^. 

The technique used to compute the physical nucleon state appearing in eq. (fj) is 
based on a double projection of the mean-field solution on linear and angular momen- 
tum eigenstates. This technique was already used to compute the static properties of 



the nucleon [16], the unpolarized and the longitudinally polarized distribution func- 
tions |l7j and the nucleon electromagnetic form factors []l9~| . We refer the reader to 
these references for more details. 

The intermediate states labelled by the quantum numbers a in eq. ([|) are also 
computed within the CCDM. Notice that they are admitted in the model, since this 
has no color. The lightest states contributing to the quark distributions are the diquark 
states (scalar ud and vectorial uu,ud). These correspond to diagrams in which a 
quark is extracted and probed by the photon. It turns out that more massive states 
(Sqlq states arising from an antiquark insertion) give smaller contributions to the 
structure functions and are important only at small x. The antiquark distribution 
receives contributions from the 4q states, which correspond to diagrams with a quark 
insertion. We explicitly found that all these terms saturate with a ~ 4% accuracy 
the normalization of the u and d valence distributions. This fulfillment of the valence 
number sum rule is of course a crucial check of the reliability of our calculation. The 
momentum sum rule is satisfied as well: in our model |14], [17]] , at Qq the valence carries 
about 75% of the energy-momentum, the remaining part being carried by the mesons 
and the dielectric field (which, in the spirit of the CCDM, embodies nonperturbative 
2 The single-minimum CCDM seems to be preferable in the light of quark matter calculations |i~8|. 
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glue). 

The transverse polarization distributions of quarks and antiquarks obtained from 
eq. (^|) using the chiral chromodielectric model are shown in Figs. 1-2. We should 
recall that the distributions computed in a quark model have no dependence on the 
momentum transfer. They represent a picture of the nucleon at some low scale Qq, 
the "model scale" . Since at such low scales higher twist effects are important, the 
structure functions obtained in quark models do not necessarily describe the physical 
nucleon at Qq, but can be used as initial conditions for the Altarelli-Parisi evolution 
from Qq to a larger scale, where higher twist contributions are absent. In previous 



works [14], [17| we showed how to determine the model scale by comparing the model 
prediction for the valence momentum with the experimental value and found for the 
CCDM Qq = 0.16 GeV 2 . We start from this scale the QCD evolution of our transverse 
polarization densities. 

Being chirally odd, hi(x,Q 2 ) does not mix with gluon distributions, which are 
chirally even. Thus its Q 2 evolution at leading order is governed only by the process 
of gluon emission. The Altarelli-Parisi equation for the QCD evolution of hi(x, Q 2 ) is 

dlog<^ 2n Jx y y 

where the leading order splitting function Ph(y) has been computed by Artru and 

Mekhfi || and reads 



P h (y) = 3 



(1 + 2/)- 



2 + £*(y-i; 



(12) 



The Mellin transforms of the splitting function Ph(y) are the anomalous dimensions 7^ 
which govern the Q 2 dependence of the moments of hi, h^i\Q 2 ) = Jq 1 dxx n ~ l h\{x, Q 2 ), 
according to the multiplicative rule 



hf\tf) = hr{Q 



(n). 



a s {Q 



(>'/ 



(«) 



a s (Q 2 



(i) 



where nj is the number of flavors. In particular, since 7^ 
of hi and the tensor charge 8q = / dx {h\ — h\) decrease with Q 2 as 



(13) 

-2/3, the first moment 



HQ 2 ) = Sq(Ql) 



a s (Q 2 ) 



-4/27 



(14) 
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Hence the Q 2 evolution of hl(x,Q 2 ) and Sq(Q 2 ) is different from that of the helicity 
distributions Aq(x, Q 2 ) and of the singlet axial charge A£(Q 2 ). The latter is constant 
in Q 2 , being related to the matrix element of a conserved current. Therefore, although 
all existing model calculations (including ours) give results for h\ very close to those 
obtained for Aq, one should keep in mind that this scenario is valid only at the model 
scale Qq. At typical experimental scales h\ and Aq are different in magnitude and shape 
as they have evolved differently from similar inputs, and the assumption h\ ~ Aq is 
no longer tenable. 

The evolved distribution functions at Q 2 = 25 GeV 2 are also shown in Figs. 1-2. 
The tensor charges at this scale are: 5u = 0.969, 5d = —0.250. 

To illustrate the different evolution of the longitudinal and the transverse polariza- 
tion distributions we compare h\ and Au in Fig. 3. It is evident that, although at Qq 
the two distributions are almost identical, after the evolution they are largely different 
at small x. In particular, the transverse distribution is considerably smaller than the 
longitudinal one for x < 0.1. A similar situation occurs for the d distributions. 

Let us turn now to the possible determination of h±. The most promising way to 
detect the transverse polarization distribution is to measure the double-spin asymmetry 
in the Drell-Yan process with two transversely polarized proton beams. This quantity 
is given by (see e.g. 

E q e 2 q hl(x aj M 2 )hi(x b ,M 2 ) + (a^b) 
TT ^ E q e 2 q q(x a ,Mi)q(x b ,M*) + (a^b) ' 1 j 

where we have labeled by a, b the two incoming protons, the virtuality M 2 of the quark 
and antiquark distributions is the squared mass of the produced dilepton pair, and 
x a ,Xb, M 2 are related to the center of mass energy yfs by x a x b = M 2 /s. The partonic 
asymmetry a^T is calculable in perturbative QCD IJ and varies between —1 and 1. 
The double longitudinal asymmetry All has an expression similar to (|15D, with the 
transverse distributions replaced by the longitudinal distributions Aq(x,M 2 ). 

In Fig. 4 we show our predictions for AttI&tt at \fs = 100 GeV 2 and for various M 2 
values. For comparison All/clll is also shown. Notice that the transverse asymmetry 
is an increasing function of the dilepton squared mass; however it remains about one 
order of magnitude smaller than the longitudinal asymmetry. In Fig. 5 we present Att 
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for x a — Xb = as a function of the center of mass energy: one can see that increasing 
y/s leads to a further depletion of A TT . The difference between A LL and A TT is an 
effect of the different evolution of h\ and Aq in the small- £ region, which dominates 
the Drell-Yan asymmetries. 

The present calculation leads us to conclude that Att is much smaller than it was 
expected on the basis of naive estimates. This is confirmed by a model-independent 
study of Drell-Yan asymmetries which will be reported in a separate paper. The 
extraction of hi is then a major challenge for experimentalists but is certainly worth 
attempting as it can add an important piece of information to our knowledge of the 
proton. 

One of us (VB) would like to thank the Institute of Nuclear Theory at the University 
of Washington for its hospitality and the U.S. Department of Energy for partial support 
during the completion of this work. 
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Figure Captions 



Fig. 1 The transverse polarization distribution of quarks hi (x) at the model scale Ql = 
0.16 GeV 2 (dashed line: h\\ dotted line: hf) and at Q 2 = 25 GeV 2 (solid line: 
h\\ dot-dashed line: hf). 

Fig. 2 Same as Fig. 1 for the antiquark distributions h\. 

Fig. 3 Comparison of the evolution of the transverse polarization distribution h\ (dashed 
line: Q 2 = Q 2 = 0.16 GeV 2 ; solid line: Q 2 = 25 GeV 2 ) and of the longitudinal 
polarization distribution Au (dotted line: Q 2 = Ql = 0.16 GeV 2 ; dot-dashed 
line: Q 2 = 25 GeV 2 ). 

Fig. 4 Predictions for the Drell-Yan double transverse asymmetry A TT /a TT (dot-dashed 
line: M 2 = 50 GeV 2 ; dashed line: M 2 = 25 GeV 2 ; solid line: M 2 = 10 GeV 2 ). For 
comparison, the double longitudinal asymmetry A^ja^ is shown for M 2 = 10 
GeV 2 (dotted line). All curves are obtained with y/s = 100 GeV. 

Fig. 5 Dependence on M 2 of the transverse double spin asymmetry at x a — Xb = (dot- 
dashed line: y^s = 100 GeV; dashed line: = 300 GeV; solid line: y/s = 500 
GeV). 
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